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If the hemispherical power asymmetry observed in the cosmic microwave background (CMB) on 
large angular scales is attributable to a superhorizon curvaton fluctuation, then the simplest model 
predicts that the primordial density fluctuations should be similarly asymmetric on all smaller scales. 
The distribution of high-redshift quasars was recently used to constrain the power asymmetry on 
scales k ~ 1.5h Mpc -1 , and the upper bound on the amplitude of the asymmetry was found 
to be a factor of six smaller than the amplitude of the asymmetry in the CMB. We show that 
it is not possible to generate an asymmetry with this scale dependence by changing the relative 
contributions of the inflaton and curvaton to the adiabatic power spectrum. Instead, we consider 
curvaton scenarios in which the curvaton decays after dark matter freezes out, thus generating 
isocurvature perturbations. If there is a superhorizon fluctuation in the curvaton field, then the rms 
amplitude of these perturbations will be asymmetric, and the asymmetry will be most apparent on 
large angular scales in the CMB. We find that it is only possible to generate the observed asymmetry 
in the CMB while satisfying the quasar constraint if the curvaton's contribution to the total dark 
matter density is small, but nonzero. The model also requires that the majority of the primordial 
power comes from fluctuations in the inflaton field. Future observations and analyses of the CMB 
will test this model because the power asymmetry generated by this model has a specific spectrum, 
and the model requires that the current upper bounds on isocurvature power are nearly saturated. 

PACS numbers: 98.80.Cq, 98.70.Vc, 98.80.Bp 



I. INTRODUCTION 

The cosmic microwave background (CMB) 0, 0, H 1, 
HUBS PJLpJIl E El 13 and the distribution 
of galaxies [l5|, Ell tell us that the early Universe was 
homogeneous on superhorizon scales, spatially flat, and 
contained a nearly scale-invariant spectrum of adiabatic 
fluctuations. These features of the early Universe pro- 
vide compelling evidence for inflation [Tt], El Ell- In- 
flation also predicts that the observable Universe should 
be statistically isotropic; any anisotropy that may have 
existed prior to inflation would be stretched beyond the 
cosmological horizon during inflation. 

There are indications, however, that the distribution 
of density perturbations is not isotrop ic J20L l2ll SI SI 
S3, SI M, S3, SI, SI, HI HI, S3, Si S3M n this arti- 
cle, we will focus our attention on one of these anoma- 
lies: the rms temperature fluctuation in the CMB on 
one side of the sky is larger than on the other side 
[H S3, SI, SE S3- This hemispherical power asymme- 
try can be parameterized as a dipolar modulation of the 
temperature anisotropy field (23 . S3]; the temperature 
fluctuation in the h direction is 



ST 

— (n) = s (n) [l + A(n-p)}, 



(1) 



where s(h) is an isotropic Gaussian random field. 1 The 
magnitude of the asymmetry is given by A and its di- 
rection is specified by p; the most recent analysis [13], 



[1] This parameterization is based on a phenomenological model pro- 
posed in Ref. [35l . 



using the WMAP5 Internal Linear Combination (ILC) 
map [H, found A = 0.072 ± 0.022 for I <, 64 with p 
pointing at (£,b) = (224°, -22°) ± 24°. No explanation 
for the asymmetry involving foregrounds or systematics 
has been forthcoming, and only a few models for a pri- 
mordial origin have been proposed [33, SI, Sil ■ 

In Ref. [39| . Erickcek, Kamionkowski, and Carroll an- 
alyzed how a superhorizon fluctuation in an inflation- 
ary field could generate such a power asymmetry. We 
found that the power asymmetry cannot be reconciled 
with single-field slow-roll inflation; the superhorizon fluc- 
tuation in the inflaton field that is required to generate 
the observed asymmetry would also induce unacceptable 
anisotropy in the CMB temperature on large angular 
scales. We then considered an alternative inflationary 
theory, the curvaton model [H S3, SI Si], which had 
been su gges ted as a possible source of a power asymme- 
try [37l |44|. In the curvaton model, the inflaton field 
dominates the Universe's energy density during inflation 
and drives the inflationary expansion, but the primordial 
fluctuations arise from quantum fluctuations in a sub- 
dominant scalar field called the curvaton. In Ref. [3^ ]. 
we showed that a superhorizon fluctuation in the cur- 
vaton field can generate the observed asymmetry while 
respecting both the homogeneity constraints imposed by 
the CMB [45[ and the constraints imposed by upper lim- 
its to non-Gaussianity [H, ES S3, HE S3 • 

The model discussed in Ref. [li] predicts that the 
magnitude and direction of the power asymmetry are 
scale- invariant. There are indications, however, that the 
asymmetry in the CMB temperature fluctuations has a 
smaller amplitude at i ~ 220 [28| and does not extend 
to I <; 600 [H S3] ■ Furthermore, an analysis of quasar 
number counts reveals that any asymmetry in the direc- 
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tion (l,b) = (225°, -27°) in the rms amplitude of pri- 
mordial density fluctuations on scales that form quasars 
(k ~ 1.3ft-1.8ft Mpc" 1 ) must correspond to A <, 0.012 at 
the 95% C.L., assuming that the perturbations are adia- 
batic [BlT ]. In this article, we consider how a superhorizon 
fluctuation in the curvaton field could produce a scale- 
dependent power asymmetry that is more pronounced on 
large scales than on small scales. 

It is possible to dilute the power asymmetry on smaller 
scales by introducing discontinuities in the inflaton po- 
tential and its derivative that change the relative con- 
tributions of the curvaton and inflaton fields to the pri- 
mordial perturbations [37t f30[ | . We examine this broken- 
scale-invariance model in Appendix [A] and find that the 
discontinuity in the inflaton potential required to sat- 
isfy the quasar constraint on the asymmetry violates con- 
straints from ringing in the power spectrum [52|, [53| . In 
Appendix [A] we also find that it is not possible to suffi- 
ciently dilute the asymmetry on small scales by smoothly 
changing the relative contributions of the curvaton and 
inflaton fluctuations to the primordial power spectrum. 

In light of these difficulties, we turn our attention to 
the dark-matter isocurvature perturbations g enerate d by 
some curvaton scenarios [H, [Hj, HE HE [H?], HE HE HE Hlf . 
In the presence of a superhorizon fluctuation in the cur- 
vaton field, the power in these isocurvature perturbations 
will be asymmetric. Since isocurvature perturbations de- 
cay once they enter the horizon, they will contribute more 
to the large-scale (I ^ 100) CMB anisotropies than to 
the smaller scales probed by quasars. Even though the 
asymmetry in the adiabatic perturbations, which is di- 
luted by the inflaton's contribution, and the asymme- 
try in the isocurvature perturbations are scale-invariant, 
the total asymmetry will be suppressed on subhorizon 
scales as the isocurvature perturbations' contribution to 
the total power decreases. Consequently, the desired 
scale-dependence of the asymmetry is a natural feature 
of isocurvature perturbations. In this article we demon- 
strate that, in certain curvaton scenarios that produce 
dark-matter isocurvature perturbations, a superhorizon 
fluctuation in the curvaton field can produce the observed 
asymmetry in the CMB without violating any other ob- 
servational constraints. 

We begin by briefly reviewing how isocurvature pertur- 
bations are generated in the curvaton scenario in Section 
HT1 and we review the CMB signatures of isocurvature 
perturbations in Section IIIII In Section IIV1 we examine 
how a hemispherical power asymmetry could be created 
by a superhorizon fluctuation in the curvaton field in two 
limiting cases of the curvaton scenario. We find in Sec- 
tion |IVX] that it is not possible to generate the observed 
hemispherical power asymmetry if the curvaton decay 
created the dark matter because the necessary super- 
horizon isocurvature fluctuation induces an unacceptably 
large temperature dipole in the CMB. In Section IIVB1 
we show that the observed asymmetry can be generated 
by a superhorizon curvaton fluctuation if the curvaton's 
contribution to the dark matter is negligible. Our model 



predicts that the asymmetry will have a specific spec- 
trum and that the current bounds on the contribution of 
isocurvature perturbations to the CMB power spectrum 
are nearly saturated. We summarize our findings and 
discuss these future tests of our model in Section [V] As 
previously discussed, we show in Appendix [A] that it is 
not possible to give the asymmetry the required scale- 
dependence by changing the relative contributions of the 
curvaton and inflaton to the adiabatic power spectrum. 
Finally, we provide a more detailed description of how 
the curvaton isocurvature perturbation can generate a 
dark-matter isocurvature fluctuation in Appendix [B] 



II. ISOCURVATURE PERTURBATIONS IN THE 
CURVATON SCENARIO 

In the curvaton scenario [IE EH, El, EH , there is a sec- 
ond scalar field present during inflation, and the energy 
density of this curvaton field is negligible compared to the 
energy density of the inflaton. The curvaton (er) is as- 
sumed to be a spectator field during inflation; it remains 
fixed at its initial value cr* and its energy is given by its 
potential V{a) = (l/2)mla 2 , with m a <C Hinf, where 
.ffinf is the Hubble parameter during inflation. When 
H ~ m a after inflation, the curvaton field begins to oscil- 
late in its potential well, and it behaves like a pressureless 
fluid until it decays. We will assume that the curvaton 
field is non-interacting prior to its decay. 

During inflation, quantum fluctuations in the curva- 
ton field [(5a) lms = H in f/(2ir)] generate a nearly scale- 
invariant spectrum of isocurvature fluctuations. After 
the inflaton decays into radiation, the growth of the cur- 
vaton energy density relative to the radiation density 
creates adiabatic perturbations from these isocurvature 
fluctuations. If the curvaton decays before any parti- 
cle species decouples from radiation, then the isocurva- 
ture fluctuation is erased after the curvaton decays be- 
cause isocurvature fluctuations between interacting fluids 
in thermal equilibrium decay quickly [EE HI] • If the cur- 
vaton decays after a particle species decouples from the 
radiation, however, there is a lasting isocurvature fluctu- 
ation between that species and the radiation in addition 
to the adiabatic perturbation generated by the growth 
of the curvaton energy densityrelative to the radiation 
density after inflation [H HI HE HE H3, HI, HE HE HI ■ 

We will restrict our attention to scenarios in which the 
curvaton decays after dark matter freeze-out but prior 
to the decoupling of any other particle species. In this 
case, an isocurvature fluctuation between dark matter 
and radiation is created. (We will neglect baryon isocur- 
vature modes, which may arise due to the annihilations of 
baryons and antibaryons created during curvaton decay 
[IEHI1])- I n this section, we will summarize how the final 
adiabatic perturbation and the dark-matter isocurvature 
fluctuation relate to the initial curvaton perturbation. A 
more detailed review of the relevant physics is presented 
in Appendix [Bl 
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Working in conformal Newtonian gauge, we take the 
perturbed Friedmann-Robertson- Walker (FRW) metric 
to be 

ds 2 = -(1 + 2*)d< 2 + a 2 (t)%(l - 2$)dx i dx j , (2) 
where a is normalized to equal one today. We define 



d = -* - H- 



Pi 



(3) 



to be the curvature perturbation on surfaces of uniform 
i-fluid density, and 



(4) 



is the curvature perturbation on surfaces of uniform to- 
tal density. Throughout this paper, a dot refers to dif- 
ferentiation with respect to proper time t. We use the 
notation S*i 7 = 3(Ci — C,^), where a subscript 7 refers 
to radiation, to describe isocurvature fluctuations. For 
any non-interacting fluid, Q is conserved on superhorizon 
scales. In the absence of isocurvature perturbations, Q is 
constant on superhorizon scales, but if there is an isocur- 
vature perturbation, then £ evolves due to the changing 
value of pi/ p. 

Immediately after inflation, there are superhorizon adi- 
abatic fluctuations from inhomogeneities in the inflaton 
field £W ~ C-y^ and superhorizon isocurvature fluctua- 
tions in the curvaton field given by S ai . After curva- 
ton decay, there are superhorizon adiabatic perturbations 
£( f ) and superhorizon dark-matter isocurvature pertur- 
bations 5 m7 . These perturbations are related through a 
transfer matrix: 



C (f) 
5, 



1117 



1~ss 



C (i) 

S, 



(77 



(5) 



This transfer matrix is completely general and applicable 
to the evolution of any mixture of isocurvature and adi- 
abatic perturbations. The left column indicates that su- 
perhorizon adiabatic perturbations do not evolve in the 
absence of isocurvature fluctuations and that they are 
incapable of generating isocurvature fluctuations. The 
expressions for T^s and Tss are model dependent. 

In the curvaton scenario, T^s depends on the fraction 
of the Universe's energy that is contained in the curvaton 
field just prior to its decay. We define 



R 



4^ 7 + 3Q CT + 3n cdr _ 



(bd) 



(6) 



where f2 7 , Vl a , and f2 c dm are the radiation energy density, 
curvaton energy density and cold-dark-matter density, 
respectively, divided by the critical density. Through- 
out this paper, quantities with a "bd" superscript are 
to be evaluated just prior to curvaton decay. We will 



assume that R <C 1 so that the curvaton never domi- 
nates the energy density of the Universe. In the limit of 
instantaneous curvaton decay with i?<l [42j , 



T 



R 
"3 



-^l bd) . 



(7) 



A numerical study of curvaton decay in the absence of 
dark matter and perturbations from the inflaton [£W = 0] 
indicates that this instant-decay expression for Tqs is ac- 
curate to within 10% provided that R is evaluated when 
H = T a /1A, where T a is the curvaton decay rate j65j . 

If the dark matter freezes out prior to curvaton decay, 
then a dark-matter isocurvature perturbation is created 
when the dark matter freezes out and when the curvaton 
decays. For J{ « 1 and instantaneous curvaton decay 



T 



ss 



(a-3)fi 



(fr) 



2(q - 2) + 

R o( bd ) 



(fr) 



(bd) 
cdm 



n 



(bd) 



b r> (bd) 



(bd) 

cdm 

R, 



(bd) 



(8) 



where quantities with an "fr" superscript are to be eval- 
uated when the dark matter freezes out. In this expres- 
sion, a = d "j m |, f . gives the dependence of the rate 
for dark matter annihilations r c d m on temperature T for 
s-wave annihilations (a ~ 21 for neutralino dark matter), 
and B m = T a -^ m /T IT is the fraction of the curvaton en- 
ergy that is turned into dark matter when the curvaton 
decays. Eq. J5]) differs slightly from the expression for 
Tss in Ref. [59(, but the two expressions are equivalent 



because Cl„ /Q 



(fr) /n ( bd ) 



n (t *> /o 

cam./ ■ 



(bd) 

cdm ' 



Numerical studies 



confirm that this expression for Tss is accurate provided 
that the decay of the curvaton does not trigger a second 
era of dark matter self-annihilation [5^ |. as discussed in 
Appendix [B] 

It will be useful to make the R dependence of Tss ex- 
plicit by defining BR = B m fl^ d ^ /f^dm to ^ e tne dark 
matter density from curvaton decay divided by the dark 
matter density prior to curvaton decay. We will also de- 



fine A = (4/3V# (bd) /ff (fr) so that 
pendix|5]). In this notation, 



\R (see Ap- 



Tss — 



(a- 3)Ai? 



2(a 



XR 



1 



1 + BR 



BR 
1 + BR. 



- R. 



(9) 

In our analysis, we will consider two limiting cases: 
BR 3> 1 (i.e. the curvaton creates nearly all the dark 
matter), and BR <c 1 (i.e. the curvaton creates an in- 
significant fraction of the dark matter). In both cases, 
we will still assume that R <C 1 . 

If the curvaton creates nearly all the dark matter so 
that BR > 1, then 



lim T S s 



(a-3)XR 
2(a - 2) + XR 



(£) + i-R. do) 
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When we recall that A < 1 is required to make the dark 
matter freeze-out prior to curvaton decay, we see that 
the first term in Eq. (|10p is proportional to X/B, which 
is much smaller than R if BR 3> 1. The first term is 
therefore negligible, and we are left with 



lim Tss = I — R- 

Bi?3>l 



(11) 



In the opposite limit, in which the curvaton's contribu- 
tion to the dark matter density is negligible, we have 



lim T S s = 



(a - 3)A 
2(a - 2) 



B - 1 



R= kR 



(12) 



The first two terms in Eq. (fT2")l are positive by definition, 
so k ^ -1. The first term is always less than 0.5 since 
A < 1, but B = (4/3)B m /fi^ n } could be much larger 
than unity since fi£im ^ !• ^ ne on ly upper limit on k 
is given by BR -C 1 which implies that k <C 1/R. 
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III. ISOCURVATURE MODES IN THE COSMIC 
MICROWAVE BACKGROUND 

Now that we have defined T^s and Tss m Eq. (JSJl, 
we can relate the early-time perturbations in the matter- 
radiation fluid to the inflaton and curvaton perturbations 
created during inflation. The power spectra of the early- 
time perturbations in the matter-radiation fluid ((^ and 
Sjorf) are the spectra that we will use as initial conditions 
to calculate the CMB power spectrum. 

Following Ref. (66[, we define 



v c (k) = ^(C (f) (fc)C (f) *(fc)>, 



V s {k) 



2^2 



C C s(k) = ^(( {t) (k)S* n (k)). 



(13) 
(14) 
(15) 



We will use a similar convention for the perturbations 
from inflation: 



A 2 [ - 

B 2 



|^<C«(fc)C (i) *(fc)), (16) 



2^2 



(S a Jk)S*(k)). (17) 



Both spectra produced during inflation are nearly flat 
(e.g. [42J), and we will assume that ~ n a ~ 1. The 
initial curvature fluctuations are created by the inflaton; 
the standard slow-roll power spectrum is 



A 2 



TT6 H 



(18) 



where e# = — Hi n f/H 2 nf is a slow-roll parameter. When 
both the radiation from inflaton decay and the curva- 
ton field are perturbed, ~ 28a*/ a*, where 5a* and 



FIG. 1: CMB power spectra for unit-amplitude initial per- 
turbations. The solid red curve is Cf: the power spectrum 
derived from "P((k) = 1. The long-dashed blue curve is C^ so : 
the power spectrum derived from Vs{k) = 1. The short- 
dashed green curve is Cf 1 : the difference between the power 
spectrum derived from V((k) = Vs{k) = C(s(k) = 1 and 

/-fad I /~*iso 



a* are evaluated at horizon exit [671 ]. For superhorizon 
fluctuations, 5a and a obey the same evolution equa- 
tion, so the ratio 5a j a is conserved (56j . Given that 
Vsa = [H in {/(2Tr)} 2 7 we have 



B 2 = 



(19) 



Since £^ is determined by the inflaton fluctuation and 
5,77 i s determined by the curvaton fluctuation, S a -y and 
are uncorrelated. From Eq. ([5]), we see that 

V c (k) = A 2 + T 2 B 2 , (20) 



Vs(k) 
C C s(k) 



r s 2 s B 2 , 

TqsTssTS ■ 



(21) 
(22) 



The CMB power spectrum may be divided into con- 
tributions from adiabatic and isocurvature perturbations 
El: 



ad 



T 2 S B 2 CT 



TqsTssB C c e 



C t = (A 2 + T 2 B 2 ) C[ 

(23) 

In this decomposition, Cf d is the CMB power spectrum 
derived from a flat spectrum of adiabatic fluctuations 
with V((k) = 1, and Cf° is the CMB power spectrum 
derived from a flat spectrum of dark-matter isocurva- 
ture perturbations with Vs(k) = 1. If both isocur- 
vature and adiabatic perturbations are present, with 
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V c (k) = Vs(k) = C cs {k) = 1, then the CMB power 
spectrum is Cf d + Cf° + Cf m . Figure Q] shows these 
three component spectra, as calculated by CMBFast 
[68| with WMAP5 best-fit cosmological parameters [l^] : 
fib = 0.0462, fi c dm = 0.233, Q A = 0.721 and H = 70.1 
km/s/Mpc. 

Figure [T] clearly shows that isocurvature perturbations 
leave a distinctive imprint on the CMB power spectrum. 
It is therefore possible to constrain the properties of 
V((k) and Vs(k) using CMB data. These constraints are 
often reported as bounds on the isocurvature fraction a 
and the correlation parameter 7: 



A 2 



T? S B 2 



T 2 B 2 

s[gn{T < sTss) A 2 + T 2 s B 2 - 



(24) 
(25) 



We will find it useful to define £ to be the fraction of 
adiabatic perturbations from the curvaton: 



T 2 S B 2 



A 2 + t 2 s b 2 ' 

with T^s = R/3. We then see that 



a 



1 + m/R 2 )Th 

7 = sign(Tss)£. 



(26) 

(27) 
(28) 



Ideally, we would like to use constraints for a and 7 
that were derived assuming only that n a d — rii so ~ 1, 
where n a( j and rii so are the spectral indices for Vq (k) and 
Vs{k) respectively. Such an analysis does not exist, al- 
though constraints have been derived for the rt a( j = rtj SO 
case and have found that n a d = n- lso ~ 1 gives the best 
fit to observations [H, [z3| - Using WMAP3 data and 
large-scale structure data, Ref. [69| found that a < 0.15 
at 95% confidence, with a slight preference for negative 
values of 7, although 7 = was included in the la in- 
terval. Ref. updated this analysis and found similar 
constraints on 7, but unfortunately they did not report 
a constraint on a. Meanwhile, the most general analyses 
[66l . l69l . l70l . l7lT [72] | make no assumptions regarding nj so 



and conclude that models with n;. 



2 — 4 provide the 



best fit to the data. Since their bounds on a and 7 are 
marginalized over rii so values that are unreachable in the 
curvaton scenario, these constraints are not applicable to 
our model. 

There are analyses that specifically target the curvaton 
scenario, but they assume that the curvaton generates 



all of the primordial fluctuations (i. e. A <C 



n s B 2 ) 

[I2> Ell UM ■ In this case, £ = 1, and the isocurvature and 
adiabatic fluctuations are completely correlated or anti- 
corrclated, depending on the sign of Tss- Furthermore, 
Eq. ([27]) shows that a £ 0.9 if £ = 1 and T S 2 S £ R 2 . Since 
this high value for a is thoroughly ruled out, these analy- 
ses of isocurvature perturbations in the curvaton scenario 



disregard the possibility that BR <C 1 and assume that 
most of the dark matter is created by curvaton decay. In 
this case, Tss is given by Eq. (TTTj) and the derived upper 
bound on a (a < 0.0041 from Ref. [l2j]) implies that 
R > 0.98. Since we require R -C 1, we can conclude that 
we will be restricted to mixed-perturbation scenarios in 
which both the curvaton and the inflaton contribute to 
the adiabatic perturbation spectrum and £ < 1. 

Finally, some analyses constrain completely uncor- 
rected (7 = 0) isocurvature and adiabatic perturba- 
tions (a.k.a. axion-type isocurvature) with nj SO = 1 
[l2l [73I [zl] . These constraints are most relevant to our 
models, however, because we will see that £ = I7I must 
be small to create an asymmetry that vanishes on small 
scales. (The discussion in the previous paragraph also 
foreshadows the fact that ^ < 1 will be necessary to ob- 
tain R <ti 1.) Furthermore, we will show that only mod- 
els with 7 > can generate the observed asymmetry, so 
the constraints derived in Refs. [69|, [7(| are too gener- 
ous. We will therefore use the bound on a derived for 
uncorrelated adiabatic and isocurvature perturbations in 
our analysis. WMAP5 data alone constrains a < 0.16 
at 95% confidence, but the upper bound on a is signif- 
icantly reduced if BAO and SN data are used to break 
a degeneracy between a and n a d- With the combined 
WMAP5+BAO+SN dataset, the 95% C.L. upper bound 
on a is [121 



a < 0.072, 



(29) 



with a best-fit value of n a d — 1. Other analyses have 
found similar bounds: a < 0.08 [z|[ and a < 0.09 [7i| at 
95% C.L. 

The other observable effect of isocurvature fluctuations 
that we must consider is non-Gaussianity [13, [6l| [ll, [7(| 
[73. Following Ref. [z|, we define through 



(30) 



where r\ is drawn from a Gaussian probability spectrum. 1 
This is analogous to the definition of /nl in terms 
of the gravitational potential for adiabatic perturba- 
tions [7a . [79l ]. Given the current upper bound on a, 

— 10 4 produces a CMB bispectrum that is equal 
in magnitude to the CMB bispectrum if /nl ^ 17 for 
purely adiabatic perturbations (76[. Furthermore, non- 
Gaussianity from isocurvature perturbations is distin- 
guishable from adiabatic non-Gaussianity through the 
scale dependence of the bispectrum [gj EH III] , and an 



[1] This definition of /jj 1 ™' differs from the definition given in Refs. 

[6ll . l75l . |77| , which define in terms of the bispectrum of 

curvature perturbations during matter domination. That 
includes information about the isocurvature fraction a and is 
consequently much smaller than the /wl defined in Eq. I|30j l 
for a given model. 
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analysis of the WMAP5 data with Minkowski functional 

found -2.5 x 10 4 < /£l 0) < 2.0 x 10 4 at 95% C.L. for 
a = 0.072 

For isocurvature perturbations from the curvaton, 



ss 



. <5ct* 



<Sct* 
ct* 



and we can set r/ = 27ss<5ct*/ct*. Thus we see that 



r (iso) 
J NL 



4Tc 



ss 



(31) 



(32) 



,(iso) 



<, 2.5xl0 4 for 



for the curvaton model. Given that | 
a = 0.072, we see that \Tss\ 10~ 5 is required if the cur- 
rent bound on isocurvature power is saturated. The cur- 
vaton also introduces non-Gaussianity in the adiabatic 
perturbations; since the fluctuations from the inflaton 
are Gaussian [801 ] . /nl for mixed perturbations from the 
inflaton and curvaton is given by (5f| [8l[ 



/nl - m 



(33) 



The current upper limit on , f NT from the CMB and large- 
scale structure is /nl & 80 [lj, [H, S3, IH, S^] . 



IV. A POWER ASYMMETRY FROM 
CURVATON ISOCURVATURE 

In an earlier article (39| . Erickcek, Kamionkowski, and 
Carroll proposed that the hemispherical power asymme- 
try in the CMB could result from a large-amplitude su- 
pcrhorizon fluctuation in the curvaton field cr, as depicted 
in Fig. O The difference between ct* on one side of the 
surface of last scatter and its average value in the ob- 
servable Universe, Act*, will introduce a power asymme- 
try AC e in the CMB through Eq. The CMB power 
spectrum depends on ct* through B 2 , as given by Eq. (fT9")) . 
and through Tqs and 755, which are functions of R. For 
R <C 1, the Universe is radiation-dominated between the 
end of inflation and the decay of the curvaton, and 



R 



\mpij v r CT 



(34) 



where m 2 ^ = G 1 is the Planck mass [42j . Differentiating 
Eq. (j2"3"|) with respect to cr* gives 



CT * 



9 1 

R 2 dT ss 
3 dR 



T 2 
2 ss 



Cf 



1T SS R- 



,dTc 



ss 



dR 



c 1 ; 



(35) 



where we have used T^s = R/3 for the curvaton scenario. 

In Ref. we assumed that the curvaton decayed 

prior to dark matter freeze-out so that no isocurvature 
perturbations are created. In this scenario, T$s = 0, and 



AC e 

~~cT 



2 Act^ 



(36) 



This power asymmetry is scale-invariant. However, if the 
curvaton also generates isocurvature perturbations, the 
power asymmetry will be scale-dependent due to the dif- 
ferences between Cf d , Cf°, and C™ 1 shown in Fig.[Tj We 
will extract this scale-dependence by defining Kg through 



AC, 



Act* 
2— Kg 



(37) 



The dipolar modulation parameter A used by Refs. [25l . 
[27| and defined in Eq. fTj describes the asymmetry in 
the amplitude of temperature fluctuations, so for small 
A, A ~ (l/2)(ACg/Cg). The modulation is assumed to 
be scale-invariant and is measured for t < £ max . To re- 
late the scale-dependent power asymmetry described by 
Kg to A, we assume that all modes with 2 < t < £ max 
are weighted equally in determining the measured asym- 
metry. Since there are (£ max — l)(^max + 3) modes in 
total, 



A = 



Act* 



21+ 1 



1=2 ^ n 



3) 



(38) 



We note that A does not depend on the amplitude of 
the superhorizon fluctuation; it is determined by T$s, R, 
and £. Since Act* cannot be larger than cr*, A is the 
largest asymmetry that can be produced by a particular 
curvaton scenario. Unless otherwise noted, we set £ max — 
64 to match Ref. [27j ■ As mentioned previously Ref. 
[13 found that A = 0.072 ± 0.022 for t <, 64, yet the 
isotropic distribution of quasars constrains A <, 0.012 for 
k ~ 1.3/j- 1.8/i Mpc" 1 [5l|. 

In the following subsections we will examine Kg for the 
two limiting cases discussed in Section HU First, we will 
consider scenarios in which most of the dark matter is 
created during curvaton decay and TgS — 1 — R. Then 
we will consider scenarios in which the curvaton's con- 
tribution to the dark matter is negligible and 755 = kR 
with —1<,k<,1/R. In both cases, we will see that 
Kg decreases rapidly when I <; 10. We will also find 
that the superhorizon curvaton fluctuation required to 
generate the observed asymmetry must have a large am- 
plitude: Act*/ct* <; 1/2. We therefore must consider 




FIG. 2: Measurements of temperature fluctuations in the 
CMB show that the rms temperature-fluctuation amplitude 
is larger in one side of the sky than in the other. We pro- 
pose that this asymmetry is generated by a large- amplitude 
fluctuation in the initial value of the curvaton field <r, . The 
fluctuation in cr, across the observable Universe is ACT,. 
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how this large-amplitude superhorizon fluctuation will 
create large-scale anisotropics in the CMB through the 
Grishchuk-Zel'dovich effect 

A superhorizon adiabatic fluctuation does not induce a 
prominent temperature dipole in the CMB due to a can- 
cellation between the intrinsic dipole and the Doppler 
dipole [ifl [83| , but this is not the case for superhorizon 
isocurvature fluctuations [H, [85j]. After matter domi- 
nation, the evolution of the potential 9 and the fluid 
velocity's dependence on are the same for adiabatic 
and isocurvature initial conditions (86j . Therefore, the 
induced Doppler dipole and the anisotropy from the in- 
tegrated Sachs- Wolfe effect will be the same for adia- 
batic and isocurvature fluctuations if the surface of last 
scatter is taken to be in the matter-dominated era. The 
only difference between the CMB dipole induced by an 
adiabatic perturbation and the CMB dipole induced by 
an isocurvature perturbation arises from the Sachs- Wolfe 
anisotropy; for adiabatic perturbations (AT/T)sw — 
^dcc/3, while (AT/T)sw = 2^dec for isocurvature per- 
turbations, where ^doc is evaluated at the time of de- 
coupling. Since we know that the integrated Sachs- Wolfe 
effect and the Doppler dipole exactly cancel the Sachs- 
Wolfe anisotropy for adiabatic perturbations, the residual 
temperature dipole for isocurvature fluctuations must be 

5*dec/3. 

If Sq is the initial matter isocurvature fluctuation set 
deep in the radiation-dominated era, then 'I'doc = —So/5 
[86[. We are considering dark-matter isocurvature fluc- 
tuations, so we have Sq — S mi fl c d m / (S^cdm + ^b), where 
S' m7 is given by Eq. (|3"Tj). We treat the superhori- 
zon fluctuation in the curvaton field as a sine wave: 
5a* = a^sin^k-x), where k <C Hq. By choosing this form 
for 5a, we have placed ourselves at the node of the sine 
wave, but the constraints we derive on 5a are not strongly 
dependent on this choice ps| . We decompose the CMB 
temperature anisotropy into multipole moments, 



(39) 



and it is not aligned with the asymmetry (87[. At least 
a portion of this anisotropy is attributable to our proper 
motion, but recent attempts to measure the peculiar ve- 
locity of the local group have found that the measured 
velocity is smaller than the velocity predicted by the 
CMB and misaligned with the temperature dipole, with 
a difference of 500 km/s [H, [8!|. We therefore take 
T> = 0.0034, which corresponds to a velocity of 500 km/s, 
to generate a conservative upper bound. 

The superhorizon fluctuation in the curvaton field will 
also generate a quadrupolar anisotropy in the CMB. The 
induced quadrupole is higher-order in Act, /a, because it 
originates from the quadratic term in S ai [see Eq. pip ]. 
Nevertheless, the upper bound on A<j*/ct* from the CMB 
quadrupole is similar to the bound from the dipole be- 
cause observations of the CMB quadrupole are not con- 
taminated by our proper motion. The CMB quadrupole 
is the sum of contributions from the superhorizon isocur- 
vature perturbation and the superhorizon adiabatic per- 
turbation (£ = T(;sS a -y during radiation domination) gen- 
erated by the curvaton field. In the coordinate system 
aligned with the asymmetry, 



»20 



(42) 



v" 1 I 



+ 5. 



ISO 

2 



cdm 



fly 



2T, 



ss 



where £f d is derived by analyzing the Sachs- Wolfe effect, 
the integrated Sachs- Wolfe effect, and the fluid velocity 
at the surface of last scatter generated by a superhorizon 
adiabatic perturbation (see Ref. [Hj]). In the limit that 
decoupling occurs after matter-domination, (5| d = 0.338, 
and S 2 so — 5/3 + <5f d since only the contribution from 
the Sachs- Wolfe effect is different for isocurvature initial 
conditions. It follows that the upper-bound on Aa»/c7» 
from the CMB quadrupole is 



and we find that, to lowest order in kxd cc , where Xd cc is 
the comoving distance to the last scattering surface, the 
dipolar moment generated by the superhorizon curvaton 
fluctuation is 



1 4lT fu \ 
aio = -r \ -T- (kx dec ) — 

O V O * 'cdm 



cdm 



ar 



(40) 



where we have chosen axes that are aligned with the 
asymmetry (z = k). The variation in a across the sur- 
face of last scattering is Aer» = a^{kxdcc), and it is con- 
strained by the CMB temperature dipole: 



Tss 



Act, 



& 0.9V, 



(41) 



where T> is the largest value of |aio| that is consistent with 
observations of the CMB dipole. The observed tempera- 
ture dipole has amplitude AT/T = (1.231 ±0.003) x 10~ 3 



(0.34 1.67 Tss) ^4.7Q, (43) 



where Q is the largest value of |a2o| that is consistent 
with observations of the CMB quadrupole. As discussed 
in Refs. [H, [45| , contributions to 020 from smaller scale 
perturbations could partially cancel the contribution to 
a2o from a superhorizon perturbation. The power in 
these fluctuations is given by the predicted value for C2 
in the best-fit ACDM model, C 2 = 1.7 x lO" 10 , and 
the WMAP5 ILC map gives a 20 = 7.3 x 10~ 6 . We 
will focus on models with T$s > 0, and we see from 
Eq. (|42[) that 020 is negative in this case. We therefore 
set Q = 1 7.3 x 10~ 6 - 2 V / C i 2~| = 1.9 x 10~ 5 as a 2a up- 
per bound on the temperature quadrupole induced by 
the variation in the curvaton field across the observable 
Universe. 
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A. Case 1: The curvaton creates most of the dark 
matter. 

If most of the dark matter is created when the curvaton 
decays, then T$s ~ 1 — R, as in Eq. (fTTj). In this case, 
Eqs. flip and {3S) imply 



2^ 



C: 



ad 



9 /Mso 



3C, C 



C; 



ad 



9 rijso 



(44) 



where we have kept only the leading-order term in R in 
the coefficients of Cf° and Cf or . We can also assume 
that RCf < Cf since Fig. Q] shows that Cf° ~ (7£ OT . 
Finally if i? £ 0.01, then i? 2 (7 £ ad < Cf for ^ ^ 1500, 
and we may neglect Cf d in the numerator. With these 
simplifications, we have ACi/Ce = —2(Aa*/a)Ke where 



K, 



94 r^iso 



94 r<iso ' 



(45) 



This approximate expression for Kt is useful because it 
only depends on £/i? 2 . It is accurate to within 1% for 
I < 1500 if R < 0.01 and accurate to within 0.1% if R < 
0.001. However, it does not have the appropriate limit 
for I -> oo; since Cf°/Cf -> in this limit, Eq. IpHjl 
shows that — ► £, but the approximate form goes to 
zero. 

Figure [3] shows the approximate form of Ki, given by 
Eq. (H]), for four values of £/i? 2 : 0.005, 0.006, 0.007, and 
0.0086. We see that Kt increases on large scales as £/R 2 
increases. On smaller scales, we see that Kt is not very 
sensitive to changes in £/i? 2 . Thus, to obtain the desired 
asymmetry on large scales and nearly no asymmetry on 
small scales, we just need to increase £/i? 2 ! Unfortu- 
nately, the upper bound on the isocurvature fraction a 
places an upper bound on £/i? 2 : 



a < 0.072 



4t < 0.0086. 
R z 



(46) 



The solid curve in Fig. [3] corresponds to £/R 2 = 0.0086 
and is therefore the maximal Kg curve that is consistent 
with the current limits on the isocurvature contribution 
to the CMB power spectrum. We also note that satisfy- 
ing the upper bound on a requires £ to be much smaller 
than R, and we have assumed that R <C 1. The adiabatic 
and isocurvature fluctuations are therefore uncorrelated. 

Figure [3] also shows that Kt peaks for £ ~ 10 and 
decreases rapidly as £ increases from 10 to 100. Further- 
more, the asymmetry nearly vanishes for larger £, so it 
will be easy to satisfy the quasar constraint. The desired 
scale-dependence comes at a cost though; the smaller val- 
ues of Kg at £ <; 20 dilute the scale-averaged asymmetry 
A. For (,/R 2 — 0.0086, the maximal asymmetry is given 
by A — 0.055. Thus we see that saturating the upper 
bound on isocurvature power (a) and setting Act* = <7* 
leads to an asymmetry that is almost Ict below the ob- 
served value. Moreover, the curvaton creates most of the 



0.1 



< 



0.05 



1 

- 


1 1 1 1 1 II 1 1 II 


1 1 1 1 1 1 


i 1 1 1 M i 

- 


- 


_^'""~~~N \ 




- 




\ \ 








\ \ 




- 


- 
























- 


\\\\ 
\\\\ 

■. \ 




- 


1 


1 i ii 


■4 

1 1 1 1 1 i 






10 



100 



1000 



FIG. 3: Kt for scenarios in which most of the dark matter 
comes from curvaton decay. The power asymmetry is given 
by ACe/Ce = —2(Aa*,/a)Ki. The solid black curve corre- 
sponds to ^/R 2 — 0.0086, which saturates the current bound 
on power from isocurvature perturbations. The lower curves 
have i/R 2 = 0.007 (long-dashed), 0.006 (short-dashed) and 
0.005 (dotted). For descending values of £,/R 2 , these curves 
correspond to asymmetry amplitudes A = 0.055, 0.045, 0.039, 
and 0.033. 



dark matter in this scenario; since Act*/ct* ~ 1 is required 
to generate sufficient asymmetry, this model requires that 
the dark matter density varies by a factor of unity across 
the observable Universe! Unsurprisingly, such a large 
isocurvature fluctuation is not consistent with the large- 
scale homogeneity of the CMB. Since T$s — 1 in this 
scenario, the CMB dipole constrains Act*/ct* ^ 3 x 10~ 3 
from Eq. ([4*1"]) . We conclude that the curvaton cannot 
generate the observed power asymmetry if the dark mat- 
ter is created during curvaton decay. 



B. Case 2: The curvaton's contribution to the dark 
matter is negligible. 

We now turn our attention to the opposite scenario, 
in which the curvaton's contribution to the dark matter 
density is insignificant. In this case, Eq. (fl"2"|) tells us that 
755 ~ kR, with —1 <, « <; 1/R. We will see, however, 
that this model can generate the observed asymmetry 
only if k <J 1.4. Thus, we will be considering scenarios in 
which the curvaton generates adiabatic and isocurvature 
fluctuations that are equal in magnitude (Tg S ~ ^5), in 
stark contrast to the scenarios considered in the previous 
section. We anticipate that generating comparable adia- 
batic and isocurvature fluctuations from the curvaton will 
be advantageous for two reasons. First, the asymmetry 
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can be partially contained in the adiabatic perturbations, 
which will make it easier to generate the observed asym- 
metry without violating the current bounds on isocurva- 
ture power. Second, the superhorizon isocurvature per- 
turbation generated by Act, will be proportional to R and 
can therefore be reduced by decreasing R. The down- 
side is that it will be difficult to make the asymmetry 
sufficiently scale-dependent to satisfy the quasar bound 
because the adiabatic perturbations are asymmetric as 
well. 

For 755 = kR, the power asymmetry generated 
by the superhorizon curvaton perturbation is given by 
ACg/Cg = 2(Aa*/a*)K e where, from Eqs. (25J and ([35]). 
we have 



G, 



ad 



9k 2 CI ! 



3kCV c 



Cf + £ (%n 2 Cf° 



3kC, c 



(47) 



We see that Kg — ► £ as £ — » 00 as expected; on small 
scales, the only source of asymmetry is the adiabatic 
power from the curvaton. We can therefore anticipate 
that the quasar constraint will place an upper bound on 
£. We also see that all the isocurvature contributions to 
the power asymmetry are proportional to k or k 2 , and 
this implies that the necessary scale-dependence of Kg 
will place a lower limit on \k\. 

Differentiating Kg with respect to £ and |k| reveals that 
increasing £ or |re| increases Kg, unless — 0.2^k^0, in 
which case the Cf° and Cj OT terms partially cancel on 
large scales, leaving Kg nearly scale-invariant. Unfortu- 
nately, the upper limit on isocurvature power places an 
upper limit on |k| and £: 



a < 0.072 n 2 i < 0.0086. 



(48) 



If we differentiate Kg with respect to |k| while keeping 
k 2 £ fixed, we find that increasing |k| decreases Kg for 
k > —0.3 and increases Kg for n < —0.3. Furthermore, 
the |k| — > 00 limit of Kg, with fixed k 2 £, is 



lim Kg 

\k\ — >oc; fixed k 2 £ 



9 (k 2 H) 6% 



Cf + 9 (k 2 C\ 



(49) 



Figure H shows Kg with k 2 £ = 0.0086 for various values 
of k. We see that as k increases from zero, the curves 
rapidly approach the dotted curve, which is Eq. (|4*9")l . 
If we could decrease k toward —00, the curves would 
approach this limit from below, but in Section|TT]we found 
that k ^ — 1, which corresponds to the lower solid curve. 

The observed asymmetry is A = 0.072 ± 0.022 for 
I & 64, which requires Kg <; 0.08 on average over this 
I range. Comparing Figs. [3] and H] reveals that it is 
much easier to generate the required asymmetry if the 
curvaton's contribution to the dark matter is insignif- 
icant because the peak in Kg is higher for 755 ~ R 
than for 7s 5 ~ 1. If we saturate the upper bound 
on isocurvature power by setting k 2 £ = 0.0086, then 
A <; 0.08 if < k < 1, and any positive value of k 
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FIG. 4: Kg for scenarios in which the curvaton's contribution 
to the dark matter density is negligible and Tss = kR. The 
power asymmetry is given by ACg/Cg = 2(A<7* /a)Kg. All of 
the curves have k 2 £ = 0.0086, which saturates the upper limit 
on isocurvature power. The top three curves have k = 0.6 
(top, solid), k = 1 (long-dashed), and k = 3 (short-dashed). 
The maximal scale-averaged asymmetries possible for these 
curves are A = 0.11, A = 0.081, and A = 0.062. The dotted 
curve is the limit as k — > 00, and it has A = 0.055. The 
bottom solid curve has « = — 1 and A = 0.043. 



has A > 0.055, which is less than la below the observed 
value if Act, /<t» ~ 1. Negative values of k are less promis- 
ing; k = — 1 maximizes Kg for negative n, and it gives 
A = 0.043 for £ max = 64. We will therefore only consider 
positive values for k for the rest of the analysis. From 
Eq. (fT2"]) we see that k can be positive only if B oc B m is 
greater than 1/2. We are therefore only considering cur- 
vaton scenarios in which the curvaton at least partially 
decays into dark matter, but its contribution to the dark 
matter density is small compared to what existed prior 
to curvaton decay. 

Figure 2] also illustrates how £ determines the small- 
scale value of Kg. As £ decreases and k 2 increases, Kg de- 
creases on small scales, and we see that k = — 1 and k = 1 
give the same small-scale value for Kg. Furthermore, the 
ratio of Kg on large scales to Kg on small scales decreases 
with increasing £. We want the asymmetry to go from 
A ~ 0.072 on large scales (£ max = 64 in the CMB) to 
A <, 0.012 on small scales (k ~ 1.3h - 1.8ft, Mpc). As 
mentioned above, the isocurvature perturbations' contri- 
bution to the total power is negligible on small scales, and 
any asymmetry is due solely to adiabatic perturbations 
from the curvaton, which implies that A = ^Aav/cr*) on 
these scales. A reduction in A from 0.072 on large scales 
to less than 0.012 on small scales therefore requires that 
A/£ <; 6. This requirement places an upper bound on £, 
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FIG. 5: The ratio A/£, with ^ max = 64, as a function of £ for 
four values of k 2 £: k 2 £ = 0.0086 (solid), 0.007 (long-dashed), 
0.006 (short-dashed), and 0.005 (dotted). Since K e -> £ as 
^ — » oo, this ratio illustrates the fractional enhancement in 
the asymmetry on large scales (£ < 64) compared to small 
scales. Since we require the asymmetry to be about 6 times 
larger on large scales than on small scales, we see that we 
require £ <5 0.016. 



as shown in Fig. [5l For k 2 £ = 0.0086, which saturates 
the upper bound on isocurvature power, we see that the 
required enhancement on large scale is attained only if 
£ ^ 0.016, and the upper bound on £ decreases with de- 
creasing a. This upper limit implies that the curvaton 
contributes only a small fraction of the adiabatic power 
(although it is a much bigger fraction than in the l~ss ~ 1 
case). The adiabatic and isocurvature fluctuations are 
therefore nearly uncorrelated. 

We now see that the required scale-dependence of the 
asymmetry limits its magnitude: for fixed a, the asym- 
metry is maximized if £ is large and k is small, but in- 
creasing £ makes the asymmetry more scale-invariant. 
Figure O summarizes the constraints on k and £. We see 
that only a limited region of the k-£ plane can produce 
asymmetries with A > 0.072 while also satisfying the up- 
per bound on isocurvature power (a < 0.072) and the 
quasar constraint (A/£ > 6). Since the quasar constraint 
is effectively an upper bound on £ and a lower bound 
on k, a tighter constraint on asymmetry in the distribu- 
tion of quasars could rule out models with A <; 0.072; 
the allowed A ^ 0.072 region is excluded if we require 
that A/t; > 16.4, which is a factor of 2.7 improvement 
on the current quasar constraint. If future observations 
reveal that A <, 0.055, however, then the allowed range 
of k values has no upper bound, and it will not be pos- 
sible to rule out these models by tightening the quasar 
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FIG. 6: The £— k parameter space for models in which the cur- 
vaton does not contribute significantly to the dark matter den- 
sity. The shaded region in the upper-right corner is excluded 
by the upper bound on isocurvature power (a < 0.072), and 
the left shaded region is excluded by the scale-invariance of 
the resulting asymmetry (A/£ < 6). The dotted curves show 
where the maximal possible asymmetry A equals the observed 
asymmetry ±lcr; the bottom shaded region cannot produce 
an asymmetry within la of the observed value. 



constraint. 

In contrast, lowering the upper bound on the isocurva- 
ture fraction a will always decrease the asymmetry am- 
plitudes that are accessible to our model. The minimum 
values of a in the A > 0.072 and A > 0.050 allowed re- 
gions are 0.054 and 0.037, respectively. If there are no 
isocurvature modes, then the Planck satellite should give 
a 95% upper limit on isocurvature power of a < 0.042, 
while an ideal, cosmic-variance-limited CMB tempera- 
ture and polarization map up to I — 2000 would con- 
strain a < 0.017 [z3|. Thus we see that Planck is capa- 
ble of ruling out our model if A <; 0.056 for I <, 64, and 
subsequent CMB measurements could test models with 
even smaller asymmetry (A <; 0.023). 

The maximum value for A is obtained when both the 
upper bound on a and the upper bound on £ are satu- 
rated; as shown in Fig. [6l A — 0.094 for £ = 0.016 and 
k = 0.74. Since A = (Aa*/a*)A, we see that A > 0.05, 
which is la below the observed value, can only be ob- 
tained if (Act*/ct„) <; 1/2. With this lower bound on 
(A(j*/<7*), the CMB dipole and quadrupole constraints 
given by Eqs. (|4"l"j) and (|43|l become upper bounds on R 
that are inversely proportional to k. With T> = 0.0034 
and Q = 1.9 x 10~ 5 , the CMB quadrupole provides a 
stronger constraint in the allowed region shown in Fig. [6] 
(k £ 0.7); for k ~ 1.4, R <, 0.00013 is required to sat- 
isfy the CMB constraints. Since Kg does not depend on 
R if 755 = kR, it will be possible to evade these con- 
straints without changing the asymmetry. (Even though 
£ depends on R, we can treat £ and R as independent 
variables because £ also depends on en and R does not.) 

The upper limit on R does have consequences for the 
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non-Gaussianity parameters, however. From Eq. (|33|) 
for /nl we see that the upper bound on R and the con- 
straint /nl & 80 implies that £ <, 0.093. Since this upper 
bound is much larger than the £ values required to gen- 
erate the necessary suppression of the power asymmetry 
on small scales, non-Gaussianity in the adiabatic per- 
turbations is not a concern. The non-Gaussianity from 
the isocurvature perturbations is bounded from below 
by the CMB dipole constraint; Eq. (|4Tj) implies an up- 
per bound on Tgg that leads directly to a lower bound 
on through Eq. (j3"2")h With the maximal variation 
in a* (Acr*/cr* = 1), the CMB dipole constraint implies 
that <; 82, independent of k. If we restrict ourselves 
to k <, 1.4, then the upper bound on R from the CMB 
quadrupole, R ^ 0.00013 for A<7*/<7* ^> 1/2, implies that 
f§£ > 1300. If a = 0.072, then these lower limits on fffi 
correspond to /nl ^ 0.1 and /nl 2.3, respectively [76| : 
both of these lower bounds on /nl are well within current 
observational limits and a re p robably beyond the reach 
of the Planck satellite [zS |9a H3l - 

Thus we see that this model for generating the power 
asymmetry does not require observable non-Gaussianity, 
unlike the purely adiabatic scenario described in Ref. 
[39l ] . We also note that significant departures from Gaus- 
sianity, should be they be observed in the future, can 
be accommodated by this model by decreasing R. This 
non-Gaussianity will include contributions from isocurva- 
ture perturbations, however, so its scale-dependence may 
be distinguishable from purely adiabatic non-Gaussianity 

ElMlzg. 



V. SUMMARY AND DISCUSSION 

A large-amplitude superhorizon fluctuation in the cur- 
vaton field can generate a hemispherical power asymme- 
try, provided that the curvaton is always a subdominant 
component of the Universe's energy density. If the cur- 
vaton decays while all particle species are still in thermal 
equilibrium with radiation, then the fluctuations in the 
curvaton field create only adiabatic perturbations, and 
the resulting asymmetry is scale- invariant (39|. Recent 
studies have revealed that the asymmetry is not scale- 
invariant; while Ad/Ce ~ 0.15 for i <, 64 27], there 
are indications that this asymmetry does not extend to 
I 600 26, 50], and an analysis of quasar number counts 
found AP(k)/P(k) <, 0.024 for k ~ 1.3ft - 1.8h Mpc" 1 . 
With the aim of explaining this scale-dependence, we 
have considered how the asymmetry produced by a large- 
amplitude curvaton fluctuation changes if the curvaton 
decays after dark matter freezes out. In this scenario, 
the curvaton produces dark-matter isocurvature pertur- 
bations in addition to adiabatic perturbations, and both 
types of perturbations have asymmetric power. Since 
isocurvature fluctuations decay after entering the hori- 
zon, their contribution to the CMB power spectrum is 
much greater on large scales (£ & 100) than on smaller 
scales, and the magnitude of the asymmetry will decrease 



accordingly. 

There are two limiting cases if the curvaton decays af- 
ter dark matter freezes out: the majority of the dark mat- 
ter can be created when the curvaton decays, or the cur- 
vaton's contribution to the dark matter density may be 
insignificant. In the first scenario, the isocurvature fluc- 
tuations from the curvaton are much larger than the adia- 
batic perturbations from the curvaton, and all of the adi- 
abatic power comes from inflaton fluctuations. Since only 
the isocurvature fluctuations are asymmetric, it is very 
difficult to generate the observed asymmetry without vi- 
olating the current bound on power from isocurvature 
modes. It is necessary to introduce an order-unity vari- 
ation in the curvaton density across the observable Uni- 
verse, and since the curvaton creates the dark matter in 
this model, this would have profound observational con- 
sequences. For instance, the resulting large-amplitude 
isocurvature perturbation induces a temperature dipole 
in the CMB that is far too large to be consistent with 
observations. We conclude that it is not possible to gen- 
erate the observed asymmetry with a superhorizon cur- 
vaton fluctuation if the curvaton creates the dark matter. 

The second scenario, in which the curvaton's contri- 
bution to the dark matter density is negligible, is far 
more promising. In this scenario, the curvaton produces 
adiabatic and isocurvature fluctuations of roughly equal 
amplitude. It is therefore slightly easier to generate the 
observed asymmetry, but the requirement that the asym- 
metry magnitude decrease by a factor of 6 between large 
and small scales limits the curvaton's contribution to the 
total adiabatic power to less than 1.6%. Consequently, 
the variation in the curvaton field across the observable 
Universe must be greater than 50% to generate the ob- 
served asymmetry. Fortunately, the amplitude of both 
the isocurvature mode and the adiabatic mode gener- 
ated by this superhorizon variation in the curvaton is 
proportional to the fraction R of the total energy density 
contained in the curvaton at the moment of its decay. We 
can therefore suppress any observational signature of the 
superhorizon curvaton fluctuation in the CMB without 
altering the asymmetry by decreasing the energy density 
of the curvaton. Decreasing R does increase the non- 
Gaussianity of the fluctuations created by the curvaton, 
but the resulting non-Gaussianity is well within the cur- 
rent observational bounds. 

We conclude that a superhorizon fluctuation in the cur- 
vaton field is capable of generating the observed asym- 
metry in the CMB while satisfying the upper bound on 
asymmetry in the quasar population if the curvaton's 
contribution to the dark matter is negligible. The curva- 
ton scenario employed by our model has two free param- 
eters: £ is the fraction of the adiabatic power that comes 
from the curvaton field, and k determines the strength 
of the isocurvature perturbation created by the curvaton 
through S m -y — kRSo-j. Both £ and n depend on other 
features of the curvaton model; £ depends on the slow- 
roll parameter e# , R, and the initial value of the curvaton 
field, while k depends on the fraction of curvaton energy 
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that is converted to dark matter and the dark matter 
density at curvaton decay. The asymmetry parameter 
A = (1/2) ACi/Ct for the WMAP5 ILC CMB map is 
0.072 ± 0.022 for £ <, 64. We find that only a narrow re- 
gion of the k-£ parameter space is capable of generating 
an asymmetry with A > 0.072, as shown in Fig.[6l so our 
model requires a fair amount of fine-tuning. If future ob- 
servations reveal that A ~ 0.05, then the allowed region 
opens up considerably and includes £ ~ with k ^> 1. 
Negative values of k cannot generate the observed asym- 
metry, however, which implies that the curvaton must at 
least partially decay into dark matter. 

The observational consequences of our model for the 
origin of the power asymmetry differ considerably from 
the predictions of the purely adiabatic model proposed 
in Ref. [39J. Whereas the purely adiabatic model re- 
quires £ <; 0.072 to generate the observed asymmetry 
and can work with £ ~ 1, our model is constrained to 
much smaller values of £ (£ ^ 0.016). Consequently, 
our model does not significantly change the tensor-scalar 
ratio [r = 16e#(l — £)] or the inflationary consistency 
relation. The dominance of the inflaton's contribution to 
the primordial fluctuations also implies that our model 
does not require detectable levels of non-Gaussianity, un- 
like the purely adiabatic model which predicts /nl <5 26 
for A — 0.072. There is one shared prediction, however; 
any primordial origin of the asymmetry predicts that the 
anisotropy in the power spectrum will produce signatures 
in the CMB polarization and temperature-polarization 
correlations [92l |. Our model also predicts a hemispheri- 
cal asymmetry in the isocurvature fraction. 

Finally, we note that this method of generating a scale- 
dependent power asymmetry through isocurvature per- 
turbations produces an asymmetry with a specific spec- 
trum. The magnitude of the resulting asymmetry peaks 
at £ ~ 10, rapidly decreases for £ = 10 — 100, and is 
nearly gone for I <; 100. Throughout this paper we 
have considered the scale-averaged asymmetry param- 
eter A with ^ max = 64; the observed value for this £ 
range in F-band is A = 0.080 ± 0.021 [l^To probe the 
scale-dependence of the asymmetry, Ref. 27J also consid- 
ered two other values of £ max in 1^-band and found that 
A = 0.119 ± 0.034 for £ max = 40 and A = 0.070 ± 0.019 
for £ maK = 80. Despite the rapid fall-off of the asymme- 
try generated by our model for £ ^ 10, it is consistent 
with these nearly scale-invariant results. For instance, 
if k = 0.75 and £ = 0.013 (a point near the middle of 
the allowed region in Fig. O, then A = 0.113(A<7*/<7*) 
for £ max = 40, A = 0.080(Aa*/a*) for £ max = 64, and 
A = 0.065(A<7*/<7„) for £ max — 80. Moreover, a search 
for asymmetry in the amplitude of the first acoustic peak 
found that A < 0.03 at the 95% C.L. for £ ~ 220 [H], 
which is consistent with the predictions of our model. 

There are also indications that the asymmetry is 
present, at least to some extent, out to £ ~ 600 [26j. Un- 
fortunately, the asymmetry parameterization employed 
in this analysis cannot be directly related to the A pa- 
rameter, so it is difficult to interpret these results. An 



analysis analogous to Refs. [25|,|27| out to higher I values 
is required to determine whether the scale-dependence of 
A predicted by our model is consistent with observations. 
Our model also predicts that at least 5.4% (3.7%) of the 
primordial power comes from isocurvature fluctuations if 
A > 0.072 (0.050) for £ <, 64. Future searches for isocur- 
vature fluctuations will therefore provide an additional 
test of our proposed origin of the CMB hemispherical 
power asymmetry; in particular, our model predicts that 
the Planck satellite should detect isocurvature perturba- 
tions if A £ 0.056 for £ <, 64 [z3- 
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APPENDIX A: A SCALE-DEPENDENT 
ASYMMETRY WITHOUT ISOCURVATURE 
MODES? 

In this appendix, we will attempt to make the power 
asymmetry generated by a superhorizon curvaton fluctu- 
ation scale-dependent without introducing isocurvature 
perturbations. As in Ref. [3^|, we define £ to be the 
fraction of the total power that comes from the curvaton 
(a): 



'PcAQ + 'PctW 



(Al) 



where <f> refers to the inflaton and £ is the curvature of 
uniform-density hypersurfaces defined by Eq. Q. The 
adiabatic power asymmetry amplitude is AP(k)/P(k) = 
2£(Aa/a). Since a has no scale dependence, the only 
way to make AP(k)/P(k) dependent on k is to make 
£ dependent on k. We need AP(k)/P(k) to go from 
2 A = 0.144 on large CMB scales (I <, 64) [23_to less 
than 0.024 on quasar scales (k ~ 1.5ft. Mpc" 1 )]^. 

We recall from Section UTT1 that V$ <j>(k) = A 2 , and 
V Cr7 (k) = (R/3) 2 B 2 , where A 2 and B 2 are defined by 
Eqs. (fTS)) and (fT9")l . Inserting these expressions into 
Eq. ifAT]) gives 



9?r 



( a * 
R 2 e H (k) \mpi 



(A2) 



ip! = G _1 . Thus we see that any scale- 
dependence in £ must originate from variation in the 
slow-roll parameter en during inflation, and we note that 
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en depends on the inflaton potential through [9S 

2 



e H (k) * e v (k) ee ^ 



It will be useful to define 



V{4>) 



v{4>) 



k=aH{4>) 



= ^(—) »' J '»f^ 
97T V TO P1/ 

because then £ has a simple form: 



e + l 



(A3) 



(A4) 



(A5) 



Even though JJ< 1 and ejy <C 1, e may take any value; 
e ^> 1 is possible if <C mpi, while e <C 1 can be obtained 
by decreasing R or e#. 

There are two ways that we can give £ the scale- 
dependence necessary to generate the observed asymme- 
try in the CMB and still satisfy the quasar constraint on 
small-scale asymmetry: we can make £ discontinuous by 
inserting a kink into V'(4>)/V(tf>), or we can choose the 
spectral indices of 7\ j(7 (fc) and V^^k) in such a way that 
£(fc) decreases smoothly as k increases. We will consider 
both approaches in this appendix. 



1. A discontinuity in £(&)? 

First, we will examine inflation models with broken 
scale-invariance. We suppose that £ is a step function: 
£ = £max for k < fc* and £ = £ m i n for fc > fc* where both 
£max and £ m ; n are constants. To achieve the necessary 
suppression of the asymmetry on small scales, we require 
that £mm/£max ^1/6 and fc* must be located somewhere 
between large CMB scales (fc ~ 0.0033/i Mpc" 1 ) and 
quasar scales (k ~ 1.5ft Mpc" 1 ) [5l|. A step function 
in £ requires a step function in e, which requires a step 
function in en- From £ m i n /£ ma x = 1/6 and Eq. (|A5[) . we 
see that 



1 



6 + 5e n 



(A6) 



Thus we see that the size of the necessary discontinuity 
in en is determined by £ max - If the curvaton dominates 
on large scales (£ ~ 1 and e > 1), then the drop in e# 
necessary to make the inflaton dominant on small scales 
is large. This drop is minimized if £ ma x is minimized, 
but £ m ax ^ 0.07 is required to generate the observed 
asymmetry [39l |. 

Now we have to worry about the shape of the power 
spectrum. Inserting a downward step in en leaves 
the curvaton perturbation spectrum unaltered, but it 
gives the inflaton perturbation spectrum an upward step. 
Therefore, we would expect that the total power on large 
scales would be smaller than the total power on small 
scales. In contrast, the primordial V((k) that fits CMB 



and large-scale-structure observations is nearly flat on 
all scales from k = 0.0001 Mpc" 1 to k = 0.2 Mpc" 1 
Furthermore, the value for erg, the fluctuation amplitude 
at 8ft Mpc -1 , derived from the CMB is consistent with 
the measurements from weak lensing observations [ll[ , so 
there can no major change in the primordial power spec- 
trum between the scales probed by the CMB and those 
probed by weak lensing. We conclude that the total pri- 
mordial power spectrum must be nearly scale-invariant. 

The total primordial power spectrum may be expressed 
in terms of £ and the curvaton power: 



P c (k) 



(A7) 



The ratio of power on large scales (fccMB < fc*) to the 
power on small scales (&q > fc*) is 



"Pc( fc CMB) 



-ff?nf( fc CMB) 



(A8) 



To compensate for the injection of inflaton power on 
scales smaller than fc* , we must introduce a discontinuity 
in V(4>) at the same <fi value as the discontinuity in V '((/>). 
Since H? n{ oc V(<f)), we see that 



^Ocmb) _ £ n 



>6, 



(A9) 



where fccMB,Q = aH- m { (</>cmb,q), is required to keep the 
primordial power spectrum scale- invariant. 

Thus we see that it is possible to hide the kick that 
the power spectrum gets when the inflaton takes over 
on small scales by introducing a drop in the inflation- 
ary energy scale that leads to a total reduction in power 
in both the curvaton and inflaton fluctuations. As the 
inflaton rolls across the discontinuity in the power spec- 
trum, the value of V(<fi) must drop by at least a factor 
of 6, and the potential on the lower side must be signifi- 
cantly flatter than the potential on the upper side; from 
e H oc (V \4>)/ 'V) 2 and Eq. (|A6]) . we conclude that 



1 



V'^q) = 1 

) 6 V 6 + 5 

'max 



(A10) 



for £min/£max = 1/6. Note that there is no lower bound 
on eh on either side of the potential break; V(<j>) can be 
as flat as we want it to be on large scales provided that it 
is even flatter on small scales. Therefore, we do not have 
to worry about constraints from the tensor-scalar ratio 
or the scalar spectral index. 

There is another concern, however. Even though the 
discontinuities in V(<j)) and V'(</>) conspire to preserve 
the flatness of the total power spectrum, the momen- 
tary interruption of slow-roll inflation that occurs as 
the inflaton field crosses the break could induce oscil- 
lations in the power spectrum localized around fc*. Ref. 
[52I ] analyzes the effects of a step in the mass to of a 
quadratic V — to 2 </> 2 /2 potential, and Ref. gener- 
alizes this analysis to steps in other inflaton potentials. 
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They restrict their analyses to breaks in the potential 
with AV/V < 0.2 because they do not want the infla- 
ton's kinetic energy to exceed its potential energy after 
the inflaton crosses the break. Even with this constraint, 
the smoothness of the observed power spectrum is very 
restrictive; CMB and large-scale-structure observations 
constrain AV/V ^ 10~ 3 at 99% confidence for fc* <, 0.1 

Mpc- 1 [13. 

The broken-scale-invariance models considered in Refs. 
[52l [53j include a small change in the perturbation am- 
plitude across the step and do not include perturbations 
from a curvaton field, but their findings are still very dis- 
couraging for our proposal. They find that the amplitude 
of the oscillations far exceeds the change in Vq (k) across 
the break, so it is reasonable to expect that any break in 
the potential will induce large oscillations, even if V^(k) 
is unaffected. Furthermore, the oscillations reach their 
maximum on scales that leave the horizon after (f> crosses 
the break. On these scales, the inflaton must dominate 
the power spectrum to suppress the asymmetry, so there 
is no hope of masking the oscillations in the inflaton spec- 
trum with the curvaton spectrum. We could hope to hide 
the oscillations induced by our model by setting fc* ^0.1 
Mpc^ 1 , but the potential change required by our model 
is so large that inflation may not resume after the in- 
flaton crosses the break. We therefore conclude that the 
discontinuity in V((f>) that is required by the broken scale- 
invariance model to satisfy the quasar constraint is highly 
unlikely to be consistent with observations. 



2. A smooth transition? 

Next we consider the possibility that £, and therefore 
ch, smoothly decrease as k increases. Given the differ- 
ence in k between large CMB scales and quasar scales, 
we would need 



trum is 



dln£ _ Aln£ 
d\nk ~ Alnfc 



L8 
6J 



= -0.29 



(All) 



if (d In £/d In k) is to be roughly constant over the scales 
of interest (0.0033/iMpc -1 <, k <, l.bh Mpc^ 1 ). 

The spectral index for £ is related to the spectral in- 
dices of the power spectrum for inflaton and curvaton 
fluctuations 42, gjl: 



dln£ _ dlnP C(T dlnT? 



din A; 



-(1-0 



dlnP c „ 



d In k d In k d in k 

= -2e H -£(-2e ff )-(l-0(-4e ff + 27? H ) 

= -(l-0(277H-2eH) (A12) 



where t]h is the slow-roll parameter 



<t>H 



(A13) 



and we note that j] H ~ (m^/ 8tt)[V" (()))/ V (()))}- e H [H. 
Meanwhile, the spectral index for the total power spec- 



n.-l S *%±3*1 , (A14) 
= -2e H -(l-Z)(2e H -2r) H ), (A15) 



and the tensor-scalar ratio is 



r = 16ejr(l-0- 



(A16) 



Therefore, the spectral index for £ depends only on £ and 
observable parameters: 



ding 
d In k 



1 + 



(A17) 



The scalar spectral index n s and the tensor-scalar ratio 
ratio r at k = 0.002 Mpc -1 are degenerate parameters; 
the blue tilt introduced by making n s > 1 may be com- 
pensated for by increasing r, which adds power on large 
scales. When n s is assumed to be scale-invariant, the 
2a upper limits on n s and r from WMAP5+BAO+SN 
are n s ^ 1.01 and r <, 0.22 [!§]■ The introduction of 
running in the scalar spectral index brings even larger 
values of n s and r into the 2er-allowed region of param- 
eter space, but the allowed range for the running index, 
a s = dn s /dlnfc, is negatively correlated with r and n s . 
This correlation makes it difficult to obtain a large neg- 
ative value for dln£/dlnfc because a s is also connected 
to the spectral index for £ through Eq. (|A15|) . 

To evaluate a s , we will need the running of the slow- 
roll parameters 



din ch 
dlnfc 

din t]h 
dlnfc 



2(e^ - T) H ), 



where 



_ ..<- P1 



_ m P1 H'(<j>)H'"(4>) 



H 2 



is a higher-order slow-roll parameter [! 
Eq. (|AT5)) . it follows that 

ding _ a s -(r/8) (vH-&/e H ) 



dlnfc 



r/4 



(A18) 
(A19) 

(A20) 
From 

(A21) 



Since £ < 1 by definition, demanding that the right-hand- 
side of Eq. (IA17|) be positive implies that the denom- 
inator on the right-hand-side of Eq. (|A21[) is positive. 
The combination of slow-roll parameters, r/H — ^h/^h, 
in the numerator could be positive or negative, but 
\r/H — £h/£h\ *C 1 is required by the slow-roll approx- 
imation. If a s is negative then the right-hand-side of 
Eq. (|A21[) can be positive only if \a s \ <§; r/8, which 
forces a s ~ 0. We conclude that a <; is required to 
make d In £/d In fc negative, and this constraint makes the 
largest allowed values for n s and r inaccessible. 
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To probe the possible evolution of £ during inflation, 
we derive a differential equation for We start with 

dhi£ 



din A; 



n s (k)-l + 2e H (k). 



(A22) 



This equation is simply Eq. (|A17[) . but with en instead 
of r. From Eq. (|A5P it follows that 



£ff (k) = e 



g(fc) 
!-£(*)' 



(A23) 



where ejv = e#(fc)/e(fc); e^v is therefore independent of 
k and can take any value. Since there are measurements 
of n s and its running, we will use these values for n s (k). 
The resulting differential equation for £(k) is 



dln£ 
dlnfc 



n s (A; )-l+Q! s In 



-2c 



N 



l-C(fe) 



(A24) 



As mentioned above, the 2a upper bound on a s from 
WMAP5+BAO+SN is negatively correlated with n s . 
We therefore consider several points, subject to the 
constraint a s > 0, on the boundary of the 2a er- 
ror ellipse in (n s ,a s ) space (see Fig. 4 in Ref. (l2l |) 
when choosing values for n s (fco = 0.002 Mpc" 1 ) and 
a s . Our aim is to minimize the ratio £min/£max = 
^(l.S/iMpc^V^O.OOSS/iMpc" 1 ), which must be less 
than 1/6 to satisfy the constraint on small-scale asym- 
metry from quasars. We find that £ m in/£max is mini- 
mized if we choose the maximal allowed value for n si 
n s {ko — 0.002 Mpc -1 ) = 1.03, which corresponds to 
a s = 0. With this parameter choice, the right-hand-side 
of Eq. (|A24[) is positive, so £ decreases monotonically as 
k increases. 

The only free parameters that remain are cn and the 
"initial" value of £ s = £(fc s ), where k s is the smallest k 
value of interest and the starting point of the numerical 
integration of Eq. (|A24j) . Together, cm and £ s determine 
eff(fc s ). From Eq. (IA23|1 . we see that e# decreases as 
decreases; since £(fc) is monotonically decreasing for 
all k, £H(k s ) is the maximum value that en will attain 
in the scale range of interest. Since we wish to maxi- 
mize — dln£/dlnfc, we want to choose the largest possi- 
ble value for cm- There are two upper bounds to consider. 
First, ejj <C 1 is required by the slow-roll approximation. 
If we insist that en < 0.5, then the maximum possible 
value for e^v is 



e N = 0.5 



1-6 



(A25) 



There is, however, an additional constraint. Since enik) 
is a monotonically decreasing function of k, the value of 
e/f(fc s > ko) sets a lower bound on the value of e#-(fco), 
provided that en is continuous. There is an upper bound 
on ejj(fco) that follows from the measured upper bound 
on the tensor-to-scalar ratio r evaluated at ko~. r(/co) < 
0.22 at 95% C.L. for a s = This upper bound 

implies that the maximum value of e^r should be 



cat 



0.014 

~1T' 



(A26) 



which is lower than Eq. (|A25[) for £ s <; 0.97. 

We integrate Eq. (|A24[) to obtain £(/c) for different 
values of £ s with €n given by both Eqs. (|A25j) and 
(|A"2"6)> . We set k s = 0.002 Mpc" 1 , and for each value 
of ejv, we find the value of £ s that gives the smallest 
value of Cmin/Cmax = f(l-5/i Mpc _1 )/^(0.0033/i Mpc" 1 ) 
with h — 0.7. We find that it is possible to obtain 
£min/£max 1/6 only if we violate the condition that 
r(feo) < 0.22. In that case, Cmi„/Cma X ^ 1/6 if 6 <, 0.25 
and eat is given by Eq. (|A25[) . It is encouraging that this 
range includes values for £ s that are greater that 0.07 
because this means that £ on large scales can be large 
enough to generate the observed asymmetry. The down- 
side is that this model has cr — 0.5 on large scales; given 
that £ s < 0.25, this means that r <; 6 on the largest scales 
for which this model applies. Since en is a monotonically 
decreasing function of k, it is not possible to make en 
small enough to satisfy r < 0.22 on the scales for which 
that bound applies (k ~ ko) and then have it smoothly 
increase to the value necessary to give sufficient variation 
in £. If we use Eq. (|A26|l to force r(k s ) < 0.22, then the 
minimal value for £ m i n /£max is only 0.56; this factor of 
2 reduction in the asymmetry between large and small 
scales is insufficient to satisfy the quasar constraint. 

In summary, it is possible for £ to be greater than 0.07 
on large CMB scales and then smoothly decrease by a 
factor of six between large CMB scales and quasar scales 
in a way that keeps the total power spectrum flat enough 
to be consistent with observations. However, the required 
values of £ and en on the largest scales are inconsistent 
with the upper bound on the tensor-to-scalar ratio on 
these scales. In order to satisfy this bound, we would 
have to discontinuously change the values of £ and en to 
suppress the tensor-scalar ratio on large scales, and we 
saw in the previous section that such discontinuities are 
problematic. 



APPENDIX B: THE DERIVATION OF T cs AND 
T ss IN THE CURVATON SCENARIO 

In this appendix, we briefly review the derivations of 
7^5 and Tss m the curvaton scenario. In the limit of 
instantaneous curvaton decay, the total curvature per- 
turbation cannot change during the decay of the curva- 
ton. We can therefore obtain T^s by equating £( f ), which 
is evaluated just after curvaton decay, to ^ bd ^ , which is 
evaluated just prior to curvaton decay. From Eq. ([4]) we 
have 



(0 



R 
~3 



4fi 7 + 3f2 CT + 30 cdn 



(bd) 



We will see below that Qdm 1S n °t ec l ua l to its initial 

(B2) 



So 



(Bl) 



value Ccdm = C^K an d we have defined 7f r through 



.(bd) 



Ai) 
'cdm 
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Since the Universe must be radiation-dominated after 
curvaton decay, the assumption that S< 1 also implies 
that fi£^ <C 1. Furthermore, we will see that 7f r <C R 
if R < 1. We thus obtain Eq. for T cs . 

If the dark matter freezes out prior to curvaton de- 
cay, then the dark-matter isocurvature perturbation is 
created in two stages. The first stage occurs at dark 
matter freeze-out |56| . The abundance of dark matter 
after freeze-out is determined by the expansion rate at 
freeze-out, so immediately after freeze-out, the hypersur- 
face of constant dark-matter density coincides with the 
hypersurface of constant total density. In the presence 
of the curvaton field, this hypersurface is not a hyper- 
surface of constant radiation density and a dark-matter 
isocurvature perturbation is created. The change in (dm 
during freeze-out is given by A£ c d m = (^fr/3)<S' cr7 where 



T fr = 



(a - 3)f£ r) 
2(a - 2) + fii fr) ' 



(B3) 



Since we assume that the Universe does not cease to be 
radiation dominated prior to curvaton decay, 




(B4) 



Between freeze-out and curvaton decay, the dark mat- 
ter is non-interacting and Ccdm is conserved; Eq. (|B2|) 



therefore relates Ccdm just prior to curvaton decay to its 
initial value. In the case that R <C 1, 7ft "C (2/3)R since 



IV -C T^IL (dark matter freezes out prior to curvaton 
decay). 

The second stage of dark-matter isocurvature creation 
occurs at curvaton decay. The change in Ccdm during 
curvaton decay is given by ^59] 



^cdn: 



r (bd) 

^cdm 



(bd) 



(bd) 



B m n 



(bd) 



o 1 



(bd) 



.(bd) 

3cdm 



(B5) 



We can relate S m7 to S ai by combining Eq. (|Bip for 

C (f) ~ and Eqs. ([Bl and dB5j) for C c ( d m ; the result 
is Eq. ||5}. This derivation of T$s neglected the possi- 
bility that the injection of dark matter particles at cur- 
vaton decay could raise the dark matter particle num- 
ber density to the point that dark matter particles begin 
to self-annihilate again. If dark matter annihilations re- 
sume after curvaton decay, then the final value of S mj 
is suppressed. It is possible to approximate the effect of 
a second stage of dark matter self-annihilation by multi- 
plying Tss by a factor of 1/(1 + T), where T = Tdm/H 
evaluated just after curvaton decay [59]. Since we are 
most interested in curvaton models in which the curva- 
ton produces very little dark matter, we will assume that 
T<1. We note, however, that our results can be easily 
adapted to cases where the self-annihilation of the dark 
matter after curvaton decay is significant. 
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